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1 Introduction 



Accelerators like ELSA, CEBAF and COSY will investigate hadron observables at a 
scale intermediate to the low-energy region where hadron phenomena mainly reflect the 
underlying principles of chiral symmetry and its dynamical breakdown, and to the large- 
momentum regime where the "strong" interaction has the appearance of being a pertur- 
bation on free-moving quarks and gluons. A theoretical description of the corresponding 
intermediate energy physics aims at an understanding of the interplay between hadronic 
degrees of freedom and their intrinsic quark substructure. At present the description 
of the intermediate energy region of several GeV, however, requires to use an effective 
parametrization of confinement within a fully relativistic formalism. Furthermore, in the 
description of hadrons as bound states of quarks nonperturbative methods are unavoid- 
able. On the other hand, for large momentum transfers the perturbative results of QCD 
should be met. In this paper we set up a further step in the development of an interpolating 
model to describe baryon structure in the intermediate energy region. 

The basic idea of this model is to parametrize complicated and/or unknown structures 
within baryons by means of two-quark correlations in the overall antisymmetrized color- 
antitriplet channels. This amounts to describing baryons effectively as bound states of 
constituent quarks and diquarks (for a recent digest of diquark models see e.g. [0]). Guided 
by the notion that a fully relativistic Faddeev equation determines baryons as bound states 
of three quarks and provides a complete and correct picture of baryons we will study the 
Bethe-Salpeter equation for a diquark-quark system. The minimal physics which has 
to be implemented in such a picture is to allow the diquark to couple to two quarks. 
Such a coupling gives rise to quark exchange between the diquark and the "third" quark 
within the baryon P|| . In this first investigation we will restrict ourselves to this minimal 
picture. Note that quark exchange is also required to reinstate the Pauli principle within 
a Faddeev approach using effective diquarks. 

A further important ingredient of our model is the effective parametrization of confine- 
ment for the constituents of the baryon, the quarks and diquarks. This is achieved by 
modifying the propagators of these particles such that no Lehmann representation can be 
found for these propagators. The price which has to be paid is an essential singularity at 
timelike infinite momentum carried by these propagators^. Note that such a singularity 
prohibits the use of dispersion relations. Nevertheless, we believe that the benefits of such 
a description, the absence of unphysical thresholds, is worth the loss of applicability of 
dispersion relations. 

Given the fact that there exists numerous models for baryons, different types of bag 
|| and soliton |||7|] models as well as non-relativisticQ [p| p~0|] and relativistic potential 



1 Recent studies of a confining interaction within the Klein-Gordon equation revealed an essen- 
tial singularity in the Isgur-Wise function of heavy quarks . 



2 For a full solution of the nonrelativistic Faddeev equation see ref. [11|. These authors, however, 
conclude: "As a result any nonrelativistic constituent quark model can at most be considered 
as a parametrization of the baryon energy levels, rather than as a dynamical model for light 
three quarks systems. Certainly it will not prove acceptable for future applications such as the 
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models ||12|| , one is tempted to argue that there is no need for another baryon model. 
A closer inspection, however, reveals that most, if not all, of these models display severe 
shortcomings when applied to intermediate energy reactions. Furthermore, at low energies 
these models describe e.g. the nucleon reasonably but not overwhelmingly well. These 
remark applies to models being so different as for example the Skyrmion versus the MIT 



bag. Hybrid models like the chiral bag |T3|,|T4[ have deepened our insight but have not been 



able to answer the question of the low-energy structure of the nucleon conclusively. A very 
recent effort to understand quark correlations in a solitonic background on the basis of 
the Nambu-Jona-Lasinio model jOJ have revealed that, at least within this model, quark 
binding energies due to a self-consistent solitonic background and due to direct two-quark 



and three-quark correlations are of the same order of magnitude ||16|| . Therefore, one might 
argue that low-energy observables are not sufficient to determine the structure of baryons 
even qualitatively. On the other hand, high precision data in a momentum regime where 
quarks certainly are dominant over collective (solitonic) effects but interactions are still 
highly non-perturbative will help to clarify issues related to the structure of baryons. The 
interpretation of these data, on the other hand, requires a not too complicated modelling 
of the corresponding reactions. In view of these remarks, we believe that there is still need 
for a baryon model adapted to the intermediate energy regime. 

This paper is organized as follows: In sect. || we describe our model: Constituent quarks 
and diquarks interact via a Yukawa interaction which represents the fact that a diquark 
may decay into two quarks. We discuss a modification of the quark and diquark propa- 
gators which allow to mimic confinement efficaciously. In sect. |3] we discuss the diquark- 
quark Bethe-Salpeter equation for the nucleon. Special emphasis is hereby put on the 
decomposition in the Dirac algebra. The components of the nucleon spinor are expanded 
in hyperspherical harmonics. Numerical solutions including scalar diquarks are presented. 
In sect. [| various nucleon form factors are calculated using Mandelstam's formalism. The 
vertex functions of external currents to quarks and diquarks are constructed from the 
corresponding Ward identities such that gauge invariance and symmetries are respected 
at the constituent level. In sect. |5| we summarize the results obtained so far and give an 
outlook. 



2 The Model 



As noted in the introduction the basic idea of our model is to parametrize complicated 
and/or unknown structures within baryons by means of two-quark correlations in the 
overall antisymmetrized color-antitriplett channels. Furthermore, it is assumed that scalar 
(0 + ) and axialvector (1 + ) diquarks are sufficient to describe baryons. As we are interested 
mainly in the nucleon we will only consider two flavors. Due to the Pauli principle the 
scalar and axialvector diquarks are then isoscalar and isovector, respectively. The minimal 
physics which has to be implemented in such a picture is to allow the diquark to couple 

description of electromagnetic form factors, hadronic decays, and other dynamical observables 
that are determined by the behaviour of the baryon wave function and are generally much 
influenced by quark-quark potential parameters." 
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to two quarks. Such a coupling gives rise to quark exchange between the diquark and 
the "third" quark within the baryon. In this first investigation we will restrict ourselves 
to this minimal picture. Note also that quark exchange is required to reinstate the Pauli 
principle on the three-quark level within a Faddeev approach using effective diquarks ||. 

Additionally, we want to modify the kinetic terms of the fields in a way that will allow 
for an effective parametrization of confinement. We will do this on the Lagrangian level 
only symbolically, and discuss the specific form of these terms in connection with the 
propagators used in the Bethe-Salpeter equation. 

The above description may be formalized with the help of the following Lagrangian: 



C = q A {x){i^d, - m q )f{-d 2 /m 2 q )q A {x) + A A {x){-d^ - m 2 s )f(-c?/m 2 s )A A (x) 



Here q(x) describes the constituent quark with mass m q . The scalar and axialvector di- 
quark field are called A and A M , respectively; their masses are denoted by m s and m a . 
The diquark field strength tensor is given by F^ u = d^A u — d u A^ + [A M , A u ] as the ax- 
ialvector diquark field is due to its isospin of non-abelian nature. Note, however, that 
we will not take into account the related self-interactions of this field. The function / 
stands as reminder that the tree-level propagators of these fields will be modified in order 
to describe their confinement. The Yukawa couplings between quarks and diquarks are 
denoted by g s and g a . Note that these interactions are renormalizable. Nevertheless, we 
will eventually substitute g s and g a by some additional momentum-dependent factors in 
order to study the influence of the fact that realistic diquarks are certainly not pointlike. 

Writing down the Yukawa interactions there are some ambiguities concerning the phases 
of the interaction term. As we have chosen a hermitian interaction term we expect the 
coupling constants to be real. The notation using g s a and g* a is nevertheless chosen in 
order to allow for a more complete discussion. 

In the Lagrangian (ffl) color indices are denoted by capital and isospin indices by small 
letters. The color coupling via the e-tensor is determined by the assumed color antitriplet 
nature of the diquark allowing for a color singlet interaction term. The matrix = r 2 is 
the antisymmetric generator of the isospin group, the matrices t% = {r 3 , 1,t x } = {rr 2 } 
are the symmetric generators. 

As we will use the Bethe-Salpeter equation in Euclidean space we will perform a Wick 
rotation thereby obtaining an Euclidean action from ([!]). We will do this by choosing a 
positive definite metric and hermitian Dirac matrices, {7^,7^} = 25^ and 7^ = 7^. 



+ 



+ 




CqUx) 



) 



(1) 
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In order to mimic confinement we will choose the function / in the Lagrangian ([]]) to be 



f(x) = l-e- d ^\ (2) 
i.e. the (Euclidean) quark and diquark propagators are given by 



S (p) = (l _ e -^ 2 +-|)/m^ (3) 

5 AB 



D AB {p) = -^—5 1 - e"^ , (4) 



Obviously, / modifies these propagators as compared to the ones for free Dirac, Klein- 
Gordon and Proca particles. It is chosen such that it removes the free particle pole at 
p 2 = —m 2 , i = q,s,a While d = 1 is the minimal way to mimic confinement, we will also 
study the cases d > 1 where the propagators are less modified for spacelike momenta. 
The propagators are free of singularities for all finite p 2 . However, they possess an 

essential singularity at timelike infinite momenta, p 2 = — oo. There exists no Lehmann 
representation for these propagators: They describe confined particles. With respect to 
our calculations the interesting property is simply the non-existence of thresholds related 
to diquark and quark production.^] 

The form of the propagators is motivated from studies in the Munzcek-Nemirovsky 



model fll7| . In this model the effective low-energy quark-quark interaction is modelled 
as (^-function in momentum space. These leads to a quark propagator without poles 
on the real p 2 axis and a behaviour for large negative p 2 similar to the one of @Q 
Recently, it has been shown that within this model diquarks are also confined if one 



includes a non-trivial quark-gluon vertex [19|. The diquark propagator does not have a 



pole and is indeed similar to the form (^,^) whereas meson masses as given by poles of 
the corresponding propagators stay almost unaltered as compared to the model with a 
tree level quark-gluon vertex. Thus, these investigations give us some confidence that an 
effective description of confinement via the parametrizations (|3],|],|5|) is indeed sensible. 
Nevertheless, for comparison we will also use in the following free tree-level propagators, 
i.e.f = l. 

As the color structure of the propagators is trivial, color indices will be suppressed in 
the following. In this paper we will assume isospin symmetry and therefore also suppress 
isospin indices. 



3 Note that these form of the propagators prohibits the use of dispersion relations. This is 
unavoidable because unitarity is violated by construction. 

4 Note that an infrared cutoff in the N JL model has similar effects |18| . 
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Fig. 1. The Bethe-Salpeter equation for the scalar and axial nucleon vertex functions. 



3 Diquark-Quark Bethe-Salpeter Equation for the Nucleon 



Using the Lagrangian ([!]) the Bethe-Salpeter equation for the diquark-quark bound state 
vertex functions is given by 



Hp) = -\9s\ 2 J ^i 5 s(- q ) l5 s(p' a )D(p' b Mp') 

/3 * f d A p' 
+ \l2 9s9a J W) 



- 4 rS(-q) l5 S(p' a )D^(p' b )^(p>), (6) 

- \9a\ 2 J 0^1 X S(-q)rS(p' a )D^(p' b )nP')- (7) 

These equations are diagrammatically represented in fig. |I| Hereby P = p a + p b denotes 
the total momentum of the bound state and p or p' is the relative momentum between 
quark and diquark: 

p = (1 - T])p a - r)p b . (8) 

The physical solution does not depend on the parameter rj e [0, 1], however, as it is known 
from previous studies of Bethe-Salpeter equations the accuracy of the numerical solution 
may depend on it (see e.g. [0). In these first exploratory calculations we will only use 
the choice which renders the equations as simple as possible, rj = 1/2. 



Using a static (momentum independent) approximation for the quark exchange, the 
Bethe-Salpeter equation (|| ^) can be reduced to a Dirac-type equation which can be 



solved almost analytically f2~I |, see also [22|. But the resulting baryon vertex functions 
are then independent of the relative momentum between quark and diquark and therefore 
not suitable for form factor calculations. 



To avoid UV divergencies in the integral equations (|6|) and (|7|) it is convenient to work with 
a finite extension of the interaction (quark exchange) in momentum space by modifying 
the exchanged quark propagator according to 

%)^%)^4^r (9) 



While it was shown in ref. [pi]] , that this choice leads to sensible results for the nucleon 



amplitudes, there is certainly the need to calculate this "form factor" , which can be related 
to the diquark Bethe-Salpeter amplitude, from a microscopic diquark model. Work in this 
direction is under progress |23[] (see also [0). 



Note that the color algebra for the eqs. (HQ) projected on the the color singlet has simply 
provided an overall factor —1. The corresponding factor for the color octet channel is |. 
As we assume the color singlet channel interaction to be attractive the color octet one is 
repulsive. Thus, using confined quarks and diquarks, no states in the octet channel can 
be generated. 

The isospin Clebsch-Gordan coefficients are given by 

and have been already absorbed in the coefficients in eqs. (|],|7|). 



Eqs. (0J^) determine the Bethe-Salpeter vertex functions $(P, p) and $ M (P, p). They are 
related to the Bethe-Salpeter wave function by multiplication with the free two-particle 
propagator, 

V(P,p) = S(P/2 + p)D(P/2-p)$(P,p), (11) 
V(P,p) = S(P/2 + p)D flu (P/2 - p)^(P,p). (12) 



3.1 Decomposition in Dirac Algebra 



As we want to describe a nucleon with positive energy we use a corresponding projection 
from the very beginning i.e. we require that the Bethe-Salpeter vertex functions and 
wave functions are proportional to a positive energy Dirac spinor for a particle with 
momentum P and helicity s, u(P, s). For later convenience we define the following Dirac- 
matrix-valued quantities p5| : 

X (P,P)= E HPp)^u(Ps), X IM (PP)= E ^(P,p)®u(P,s), (13) 

s=±l/2 s=±l/2 
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uj(P,p)= J2 y(P,p)®u(P,s), oj»{P,p) = W*(P,p)®u(P,s). (14) 

s=±l/2 5=±l/2 

These are then proportional to the positive energy projector 

A+ = Yl u{P,s)®u{P,s) = hl + f) (15) 

s=±l/2 

where P M = P^/iM is the unit momentum vector parallel to the bound state momentum. 
Using the projector property (A + ) 2 = A + the vertex functions can be written as 

$(P,p) = X(P,P)A + «(P, S ), (16) 
^(P,p) = X fl (P,p)A + u(P,s). (17) 



The most general Dirac decomposition for the scalar function x(PjP) involves four inde- 
pendent Lorentz scalar functions, e.g. 

X (P,p) = a x {P,p) + a 2 (P,p)P + a 3 (P,p)^> + a 4 {P,p)±{fl> -ffl), (18) 

whereas the one for x M (P,p) decomposes into twelve Lorentz axialvector functions, e.g. 

1^{P,P) = h(P,p)P» + b 2 (P,p)p» + b 3 (P,p)Y + b 4 (P,p)P»P + b 5 (P,p)p^ +... (19) 

Projection onto positive energies leads to two and six independent components, respec- 
tively. It is advantageous to choose the Dirac matrices in these expansions to be eigen- 
functions to A + . In addition, the expansion coefficients are required to fulfill the Dirac 
equation for a particle of mass ±zM. This leads to 

X {P,p) = S 1 {P,p)A+ + S 2 (P,p)EA + (20) 

X»(P,p) = A 1 (P,p) l5 P^A + + A 2 (P,p) l5 P^A+ 
+ B 1 (P,p) 7B p£SA+ + P 2 (P,p) 75 ^A + 

+ c^pph^Pn - Y) - p£2)a+ + c 2 (p,p) l5 (t(Pn + y)s - ^)a + 

(21) 

where 



E=-iy-P-p)/y/j%& (22) 

is a Dirac matrix whose interpretation is obvious in the center-of-mass frame of the bound 
state, see below, and 



p»=p»-pv(p. p ) and p£ = p£/Vl#l# ( 23 ) 

are projections of the relative momentum on the direction transverse to the total mo- 
mentum, pj,P^ = 0. Note that the normalization applied in eq. (^) leads to easily in- 
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terpretable expressions of the amplitudes in the rest frame. Nevertheless, for numerical 
treatment, we will modify the normalization slightly, see below. 



The notation in eqs. (|20|), fl2lD is such that an index 1 refers to an operator in the expansion 



being a solution to a Dirac equation with mass +%M (upper component in the rest frame 
of the bound state) and an index 2 refers to an operator in the expansion being a solution 
to a Dirac equation with mass —iM (lower component in the rest frame of the bound 
state). 

Let us detail this for the case of the scalar diquark. Obviously, in the rest frame of the 
bound state, P = (0, 0, 0, iM), the operators read 

A + =fl ?V HAW" °V (24) 



o oy ypa- o 

Thus, multipliying the amplitude \ with an positive energy spinor u, using the definition 
of x (Pf) and the relation u(P, s)u(P,t) = jj5 st = 5 st (the last equal sign is only valid 
because we have specialized to the rest frame of the bound state) one obtains: 

•wM^&ri)- (26) 

As anticipated this is the spinor for a | + particle. The more complicated case of the vector 
spinor $ M relating to the amplitude with the axialvector diquark as a constituent is given 
in Appendix A. 



3.2 Expansion in Hyperspherical Coordinates 



In order to obtain a numerical solution of the diquark-quark Bethe-Salpeter equation 
we expand the amplitudes Si,Ai,P>i and Cj, i — 1,2, in hyperspherical harmonics. The 
basic idea hereby is the would-be 0(4) symmetry of the Bethe-Salpeter equation if the 
exchanged particle would be massless. Thus, the hope is that only a few orders are suf- 
ficient for a numerically precise solution which will be indeed the case, see [2|j] and the 
next section. 



First, let us recall that in four dimensional spherical coordinates, 

p^ = p(cos <f> sin 6 sin if> , sin <p sin 6 sin x/j, cos # sin-?/;, cos if}), 
the spherical spinor harmonics are given by 



Z n jl m {lfj,6,(f)) 



\ 



7T(n + 1 + 1)1 



(26) 



Hereby the C^t^cosip) are the Gegenbauer polynomials [P^] . The ones with 0(3) angular 
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momentum / = are identical to the Chebychev polynomials of the second kind, 
n \ ( i\ rp ( I \ sin(n+l)V> 

Cft(cosV') = T n (cos^) = , (27) 

obeying the orthogonality relation 

TV 

Jdip sin 2 ^ T n (cos?/')T m (cos^) = ^5 nm (28) 
o 

The 0(3) spherical spinor harmonics yi m {9, <P) are the usual ones familiar from relativistic 
quantum mechanics. They obey the relation 

y J m/ 2,Jp) = -p<ryhi/2,m(p)- (29) 

Using the usual spherical harmonics Y/ m the explicit expression is given by 



1 / ±V^ ± m + hYi m _l/2 I 

y m (o, 0) = -_= I T " l ' m 1/2 (30) 



for j = I + | and j = I — |, / > 0, respectively. 



The amplitudes in eqs. (]20|) , (pi|) Si,A i: Bi and Cj, i = 1,2, are expanded in Chebychev 
polynomials, e.g. 

oo 

= E^.(P 2 ,p 2 )T„(P-p) (31) 

n=0 

and analogous for the functions A4, Bi and Cj. 



5. 5 Numerical method 



For the numerical solution of the Bethe-Salpeter equation we will restrict ourselves in 
the following to + diquarks, i.e the solution of eq. @ with g a = 0. The inclusion of 
1 + diquarks and accordingly the axialvector amplitudes of the nucleon is much more 
involved and will be treated in a forthcoming publication. Furthermore we choose equal 
masses for the quark and the diquark, m q = m s . Then all dimensional quantities in the 
Bethe-Salpeter equation can be expressed in units of the constituent mass m q . To allow 



a comparison of our results with the ones reported in [|24J] we also choose, if not stated 
otherwise, A = 2m q , which basically fixes the width of the interaction in momentum 
space. 

When working in the rest frame of the bound state we parametrize the relative momenta 
in eq. (||) according to 

p' = p'(cos (f) 1 sin 6' sin if)' , sin cj) sin 9' sin ip' , cos 9' sin ip 1 , cos ip'), (32) 
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Pn = p(0, 0, smip, cos^). 



(33) 



By changing the normalization for the Dirac matrix in (22) to 



5 = - P ■ p) /y/pV (34) 

and denoting the spatial part of the normalized 4— vector p by p, the generic structure of 
the integral equation can be written as 

/ is^p) o\_ _, rW_ , ( xs x {p,p>) o\ 

{smiPa 3 S 2 (P,p) 0j~ 9s J (2ttY {P,Pj n \p'aS 2 (P,p') 0j> [M) 
K{p,p',P) ~ l5 S(-q) l5 S(p' a )D(p' b ) (36) 

Using the expansion of the amplitudes in terms of Gegenbauer polynomials (|31]) and 
applying the orthogonality relation (^3), we extract the expansion functions (amplitudes) 
Si,m{p), (i — 1, 2; m — 0..oo) and obtain after integration over 0' 



s 2m (p) 



By taking the spin degeneracy into account we reduced the problem to two coupled integral 
equations, where Kij, = 1, 2) now denote real scalar quantities and can be calculated 
straightforwardly. The remaining angular integrals are given by 

7T 7T 7T 

JdY> = Jdip sin 2 ?/; [dip sinV Jd& sin 6'. (38) 
00 

For the actual solution the Gegenbauer expansion has to be terminated at a finite n = 
m max and due to the structure of eq. ( |37D m is in the range < m < m max . Equation 
P7| ) is now the most convenient expression of the Bethe-Salpeter equation for numerical 
treatment. While the B' integration can be done analytically, we perform the ip and 
ip' integration numerically by means of Gaussian quadratures with appropriate weights 
(w(z) = a/1 — z 2 , z = cos ip). The final step is then to discretize p and p' with k max grid 
points to generate a k max x k max mesh in momentum space. Actually we map the interval 
[0, 00] to [—1, 1] and determine the grid points according to a Gaussian quadrature with 
Gauss-Legendre weights. In this way a (2-m max -k max x 2-m max -k max ) matrix is generated 
which can be treated as an eigenvalue problem; for a given bound state mass M (which 
appears nonlinearly in the integral equation) we seek the corresponding coupling constant 
g s . In order to test the numerics we solved the eigenvalue problem by straightforward 
diagonalization as well as by an iteration method and found, within numerical accuracy, 
identical results. When applying the iteration method one starts with initial guesses for 
the amplitudes and iterates the matrix equation until g s converges. 

As a further check of our numerical procedure we also solved the Bethe-Salpeter equation 
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involving explicitly the wavefunction ^(P,p) by iteration, 

*(P,p) = -g 2 s J -^H(p,p>)y(P,p>) (39) 
*(P,p) = D(P/2 - p)S(P/2 + p)${P,p). (40) 

Note that the corresponding kernel H ~ 75>S'(— (7)75 is independent of P, or M, which 
in principle allows a fast determination of M with given g s . First we use an adequate 
decomposition for *f?(P,p) as in eq. (p5|) , 

• (Pp) = C 1Gl(P ' P) 1 (41) 

We then also expand the free two-particle propagator D(P/2 — p)S(P/2 + p) in terms 
of Gegenbauer polynomials (obtained in the case of tree-level propagators analytically, in 
the case of confining propagators numerically). This defines the coefficients Uij via the 
relation 

D(Pv)S(Pv)( 15l(P ' P) VfW U " u U ™ u 

T ^*iw$)- (42) 

The final step is the projection onto the expansion coefficients S im (p) and G im (p),i = 1,2, 
to yield 



^2m(p), 



t)-(tS) (43) 



G ln (p)\ _ m ^ m+ ^, kl+m - n fu k \ v&\ (s lm (p)\ (M) 

G2n( P )) to to W U k 2 \)-\S 2m {p)) 

\k\ — m\ < n < ki + m 
k\ + m = n mod 2, 

which is suitable for an iteration procedure. Note, that we used the following addition 
theorem for Gegenbauer polynomials: 

T n (z)T m (z) = T\ n _ m \(z) + T| n _ m | +2 (z) + • • • + T n+m (z). (45) 



Once the kernel H(p,p') was computed on a momentum mesh we iterated $(P, p) (ex- 
panded up to order m max ) taking into account \l/(P,p) up to order n max , as expressed in 

In the limit n max — > 00 both methods are by construction equivalent. In the actual 
calculations n max = m max + 3 proved to be sufficient as the eigenvalues for a given m max 
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quickly converged due to the asymptotic behaviour of the kernel for higher orders m and 
n. 

Before discussing our numerical results we want to mention that the linear Bethe-Salpeter 
equation provides no a priori normalization condition. While a physical normalization of 
the bound state amplitudes will be discussed in the next section here we simply demand 
5*10 (pi) — 1 (Pi denotes the first grid point of the momentum mesh), which provides a 
preliminary overall normalization. 



3.4 Numerical results 



In this subsection we discuss our numerical results. By examining the convergence prop- 
erties of the integral equation we observed that for k max > 30 the eigenvalues and eigen- 
vectors do not depend any more on the size of the momentum mesh. The numerically 
obtained amplitudes are shown in figure ^. Type I denotes the calculations where quark 
and diquark propagators with poles are included (cf. f = 1 in eq. (||)); type II denotes 
calculations with confining propagators. When doing the calculation with propagators of 
type II we also vary d, the damping factor in the exponentials. As it is seen in the figures 
if d is increased, the amplitudes approach the results obtained with type I propagators. 
For the case d = 10 there is no visible deviation from the type I amplitudes. This feature 
is expected since for a large damping factor basically only the timelike properties of the 
quark and diquark propagators are affected by the exponential. We note that the ampli- 



tudes for the type I calculations are in agreement with the results in ||24|| , which provides 
a check of our calculation. 



type I 


TTlmax 


M = m q 


M = 1.9m q 


M = 1.99m q 





15.8164 


9.9237 


8.2132 


1 


16.1141 


9.8853 


8.1729 


2 


16.0809 


9.8488 


8.1496 


3 


16.0812 


9.8485 


8.1494 



Table 1 

Eigenvalues of the Bethe-Salpeter equation with type I propagators for different nucleon masses 
when the expansion in terms of Gegenbauer polynomials is terminated at m — rn ma x. 



From figure (0) it can now be seen that the employed expansion in terms of Gegenbauer 
polynomials converge rapidly: The absolute value of subsequent amplitudes (compare e.g. 
Sn with S10) decreases roughly by one order of magnitude. So higher orders {n > 3) can 
be safely neglected. 
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Fig. 2. Bethe-Salpeter amplitudes Sio(p)...S22(p) calculated with type I and type II propagators 
(with different values for d) and with k max = 30. 
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type II, d=l 




M = m a 


M = 1.9m„ 

y 


M = 1.99m„ 

y 





17.7627 


14.0058 


13.5747 


1 


18.2856 


14.0938 


13.6373 


2 


18.2502 


14.0485 


13.5916 


3 


18.2507 


14.0487 


13.5918 



Table 2 

Eigenvalues of the Bethe-Salpeter equation with type II and d = 1 propagators for different 
bound state masses when the expansion in terms of Gegenbauer polynomials is terminated at 

Ul — ITlmax' 

In order to stress this fact, we show in table ([!]) and (0) how the eigenvalues of the 
integral equation converge when more orders of the Gegenbauer expansion are included 
in the calculation, i. e. m max is increased. The same behaviour can be observed when the 
eigenvectors (amplitudes) are examined; graphically they then become undistinguishable. 
Thus, we conclude that especially in the "physical region" of the binding energy (m q < 
M < 1.9m q ) sufficient convergence is achieved with n max = 2, thereby including 3 orders 
in the expansion. 

Though the two methods to solve the Bethe-Salpeter equation (one for the nucleon vertex 
function and one for the wave function) mentioned in the last section are expected to 
give equivalent results, deviations might arise in the numerical integration procedure. 
The kernel K(p,p'; P), including the free two-particle propagator approaches zero much 
faster than H(p,p') for large values of p' and for higher orders m,n in the Gegenbauer 
expansion. Therefore, using the second method, the number of grid points has to be 
increased to obtain the same numerical accuracy; we typically chose 50 points for the 
momentum mesh. Then the results agree very well within numerical accuracy. 

Having in mind that S% (p) describes the lower component of the nucleon Dirac field (which 
is usually neglected, e.g. in nonrelativistic calculations), it is surprising that we obtain 
rather large "small" components. This was also observed in ref. [24]. To study this in 
greater detail we plot in figure (§) the amplitude S20 (p) for different values of the bound 
state mass. For both types of propagators the maximum of the amplitudes 5*20 {p) increases 
when the binding energy Eb = 2m q — M tends to 0. Note however that this result depends 
on the chosen representation of the fermion field (in our case the Dirac representation )Q. 
We will investigate in the next section whether the large lower component of the nucleon 
spinor has effects on observables where a 75 is involved in the coupling to the external 
current (<m(Q 2 ) and g n NN(Q 2 )) because this 7-matrix induces a mixing between upper 
and lower components which give rise to contributions to these observables. 



The dependence of the amplitudes on the width of the interaction can be seen from fig. 
(|4]). Although we restricted ourselves in the plot to a calculation of 5io(p) for type II 

5 We thank F. Lenz for this remark. 



15 




Fig. 3. The lower component S2o(p) for M = m q ,1.9m q ,1.99m q using type I propagators, left 
hand figure, and type II propagators with d = 1, right hand figure. 

S 10 (P) 




0.0 2.0 4.0 6.0 8.0 10.0 



Fig. 4. S\o(p) is plotted for type II and d = 1 propagators for different values of A 2 (in units of 
m q ) which controls the range of the interaction. 

propagators with d = 1 and a bound state mass of M = 1.9m q the observed behaviour is 
quite generic: When A is increased, higher modes of the interaction are included, leading 
to a broader amplitude in momentum space. Correspondingly, in coordinate space the 
interaction become shorter ranged leading to a narrower amplitude. Although this be- 
haviour would allow to fit A e.g. to the mean square radius of the nucleon, we will not do 
this in our following calculation. 

In order to demonstrate the main advantage of our approach, namely working not only 
with free (type I) propagators but also with confining propagators (type II), we display 
in figure || the calculated coupling constant g s (~ eigenvalue) as a function of the nucleon 
mass. For type I propagators the eigenvalues near threshold start to decrease rapidly. 
Ideally they should display an infinite slope for M = 2m q , because there the bound 
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Eigenvalue vs. Bound State Mass 



16.0 
14.0 
12.0 
10.0 
8.0 
6.0 
4.0 

1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 

M/m q 

Fig. 5. Eigenvalues of the integral equation as functions of the bound state mass M. In case of 
type I propagators, the influence of the quark-diquark threshold is clearly visible. 

state decays into its constituents. To observe this feature a much more refined numerical 
method (suitable for extremely loose bound states) would be necessary. Our numerical 
method, which is very robust in a wide range of bound state masses cannot account for 
this property, but signals the nearby threshold by a strong sensitivity of g s on M. When 
using now the confining propagators we get a totally different behaviour: The eigenvalues 
are not affected by the "threshold" any more (which in this case is no threshold at all). 
The function g s is a smooth function and nothing signals the possibility of the decay 
of the bound state. In this way confinement is realized in our approach: By choosing 
confining propagators for the constituents of the nucleon, we exclude its unphysical decay 
into free quarks and diquarks. Even when the damping factor d is increased, so that for 
M < 1.9 the eigenvalues and eigenvectors of a type I calculation are basically recovered, 
near "threshold" we again observe confinement. 

Before we describe in the next section the application of the so far obtained amplitudes 
in a calculation of baryonic matrix elements, we want to mention that the numerical 
solutions of the Bethe-Salpeter equation can be parametrized for further applications 
very effectively. By examining the infrared and ultraviolet asymptotics of the amplitudes 
we find simple rational functions with parameters fitted to the numerical solutions. For 
details we refer to appendix |B[ 



4 Nucleon Form Factors 

In order to test the validity of the proposed quark-diquark picture of the nucleon, it is 
not enough just to solve the corresponding Bethe-Salpeter equation, since its solution, 
the Bethe-Salpeter amplitudes (or the vertex functions) have no physical interpretation 
by themselves. It is therefore necessary to use the solution of the Bethe-Salpeter equa- 
tion to calculate physical observables as, e.g. form factors. Especially we are interested in 
electromagnetic form factors of proton and neutron, the pion-nucleon form factor g^NN 



threshold • 




type I 

type II, d=1 
type II, d=10 
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and the axial form factor g^. These observables are accessible by coupling an appropriate 
probing external current to the nucleon. 

In ref. |]28|j29[1 a comparable calculation to our approach was performed in the context of 
the NJL-model (with unconfmed quarks and diquarks), but only static properties (vanish- 
ing momentum transfer) of the nucleon were considered; in ref. |3(| spacelike form factors 
were calulated (also within the NJL model) but assuming a static quark exchange be- 
tween quark and diquark and therefore working with a "trivial" nucleon vertex function. 
A calculation within the Salpeter approach, instead of the here employed fully relativistic 
Bethe-Salpter equation, is reported in ref. [31] and [32]. So far, there exist no work in 



which a) the BS equation was solved without any approximation for the quark exchange 
and b) the resulting solution was used to evaluate observables for finite momentum trans- 

2M for a calculation of nucleon structure functions 



however, ref. 



fer to the nucleon. See 
within this approach. 

At this point we make the following remarks: The so far developed picture of a nucleon 
consisting of a quark and a scalar diquark with equal masses is certainly to naive for a 
correct descripton of the nucleon and therefore quantitatively agreement with experimen- 
tal results cannot be expected. Nevertheless we believe that the reported studies are a 
necessary step in the development of a complete baryon model. Therefore we not even try 
to fit e.g. A or d to experimental numbers. Our aim is rather to observe, whether certain 
nucleon properties can be described reasonably. 

The calculation of hadronic matrix elements between bound state vertex functions is con- 
veniently performed in Mandelstam's formalism [33], which will be introduced in the next 
subsection. 



4-1 Mandelstam's Formalism 



In order to calculate the matrix element of an external current operator (generically 
denoted by O) between bound state amplitudes we use Mandelstam's formalism [|3j| (see 
also ]p9fl), which reads in momentum space 

(6 M ) = (NiP^Sf^O^NiPuSt)) = 1 Joo(P/,Pi), (46) 

Here Pi,Pi and Pj = Pi + Q, pf denote the total and the relative momenta of the quark- 
diquark bound state before and after the interaction with the external current, while ^ 
and ^ are the properly normalized Bethe-Salpeter wave functions. By writing explicitly 
the quark and diquark propagators, which are included in the definition of the BS-wave 
function, we obtain Mandelstam's prescription in a representation, which contains the 
nucleon vertex function, calculated in the previous section. 

Wf,Pi)= J^f^x(Pf,Pf)D(kPf-Pf)S(?Pf+Pf)x 

T 6 (p f , P f ;p h Pi)D{\P - Pi )S{\P + pMP, Pi). (48) 
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Fig. 6. The two diagrams which contribute to the nucleon matrix elements in impulse approxi- 
mation. 

This representation allows the direct evaluation of the nucleon current. The 5-point func- 
tion Tq appearing in eq. fl4"T|) and (^8|) describes the coupling of the probing current with 
momentum Q to the constituents of the bound state. 

In the following we work in a generalized impulse approximation, where the coupling 
to the diquark and the constituent quark is considered, while we neglect, in this first 
calculation, the coupling to the exchanged quark. 
Accordingly the 5-point function is defined by 



r 6 {p f , P f ; pi, Pi) = 5(p f - Pi - \Q)D-\\P t - Pi)Tl(p f , P f ; p u P) 
+S(p f -Pi+ \Q)S'\\P + Pi)V%{p f , P f ; Pi , P). 



(49) 



Inserting Tq into eq. fl48|), one observes that the matrix element is split into a quark and 
a diquark part 



(50) 



If the external nucleon legs are on their mass shells (Pf = Pf = —M 2 ), the nucleon 
current depends only on the squared momentum Q 2 of the probing external current. The 
two terms appearing in flSOl) are given by the loop integrals 



and 



JUQ 2 ) = /7^x( J P/,P/)5b + )r^ ) (p + ,p_)5(p_) J D(^) x (p l , Pl 



(27T) 



(51) 



■JUQ 2 ) = J 



(Jip x(Pf,P f )D{p + )TUp + ,p_)D(p_)S(p q ) X (Pi,Pi)- 



(2tt) 



(52) 



In appendix ^ we discuss the momentum rooting in the loop integrals and their numerical 
evaluation. 

Since the Bethe-Salpeter equation does not provide the overall normalization of the bound 
state vertex function, a physical normalization condition has to be imposed. The usual 
requirement, a unit residue of the bound state propagator at the mass pole pi, which in 
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our formalism reads 



fj^mP)[ W SGP + p)D$P-p)] X (P,p)=2A+Pi>, (53) 

then leads to the correct normalization. Since the quark exchange is independent of the 
cm. momentum P, it does not contribute to the normalization of the vertex function. 
Note that the adjoint vertex function x(P, p) has to fulfill 

X(P P ) = A + X(PP) (54) 
and is therefore given by x(P, p) = (P, p)l4- 



4-2 Electromagnetic Form Factors 



To determine the electromagnetic (e.m.) form factors one has to calculate the e.m. nu- 
cleon current. The probing current is an external photon and therefore the quark-photon 
and the diquark-photon vertex functions are needed to apply Mandelstam's formalism. 
Electromagnetic gauge invariance is manifest in the Ward-Takahashi identity, 



QJV(p + Q,p) = S~ l (p + Q)- S~ L (p), 



(55) 



which provides the connection between the longitudinal part of the (di-) quark-photon 
vertexfunction and the inverse (di-) quark propagator, and for the limit Q 2 — > in the 
Ward identity 



r M (p,p) 



_d_ 

dp. 



s-\ P ). 



(56) 



Using vertex functions in accordance with these identities means to respect gauge invari- 
ance at the constituent level. 

A vertex function which solves the above identites for the quark-photon coupling is the 
Ball-Chiu vertex IB3J3H1, 



(r«)e.m.(p,*) 



+ % 



(p + jOg 

p2 _ fc2 



A(p 2 ) + A{k 2 



(A(p 2 ) - A{k 2 )) 



"7m 

7P + yk 



- l (B(p 2 )-B(k 2 )) 



(57) 



which is determined by the quark self-energy functions A(jp 2 ) and B(p 2 ) appearing in the 
quark propagator 



S(p) = -iTpa v (j?) + a s (p 2 ) = [ij-pA(p 2 ) + B(p 2 )}- 1 



(5* 



and can be read off from eq. (|3|). The vertex function (|57| ) not only satisfies eqs. (p^) and 
(|56), but is also free of kinematical singularities, has the same properties under P,C,T 
transformations as the perturbative vertex iy^ and reduces itself to this vertex in case of 
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bare quarks (this is the case when a type I quark propagator is used). In the following we 
will restrict ourselves to the pure longitudinal Ball-Chiu vertex and will not consider any 
transversal part, which cannot be determined using e.m. gauge invariance^]. 
The vertex function of a nonperturbative (extended) scalar particle, also has to be chosen 
to fulfill eqs. fl5"5|) and (|6|); this is discussed in ref. PBfl . Again we will use the simplest 
longitudinal diquark-photon vertexfunction 



(rj)e.m.(p,*) = -(p + *), 



p 2 C{p 2 ) - k 2 C(k 2 ) C(p 2 ) - C{k 2 ) 



p 2 — k 2 p 2 — k 2 

which is determined by the self-energy function C(p 2 ) of the diquark propagator 



(59) 



Dip) = -4t^~2 = - ((P 2 + ™')CV)r • (60) 
p z + m s v J 

The explicit expression can be easily obtained when comparing with eq. In particular, 
any transversal part of the vertex function, which would describe an anomalous magnetic 
moment of the scalar diquark will not be considered. 
To calculate the electromagnetic form factors of the nucleon we insert 

r^(p + , P _) = g ? (r^) e . m .(p + ,p_) (6i) 

and 

rJ(p + ,p_) = Q < ,(rJ)e. m .(p + ,p_) (62) 

with the charge matrices in isospin space given by 

Q 9 = ±(§l + r 2 ), Q d = \. (63) 



into eqs. (|5T|) and (0), respectively, and evaluate the loop integrals. As it is well known, 
using Lorentz invariance, invariance under P,C,T and the fact that the incoming and 
outgoing nucleons are onshell, the longitudinal e.m. nucleon current can be decomposed 
into 



j;- m iQ 2 ) = A + (P /; S f )[iM B (F e (Q 2 ) - F m (Q 2 ))^ 2 + F m (Q 2 ) 7/ JA + (P,, Si). 



(64) 



Note that within our formalism the Dirac spinors usually appearing on the right hand 
side of this equation are replaced by the A + projectors. The Lorentz invariant functions 
F e (Q 2 ) and F m (Q 2 ) denote the electric and magnetic form factor, which can be extracted 
by taking appropriate traces. Here and in the following subsections the calculations are 
most conveniently performed in the Breit-frame (see appendix |C|) . 



In ref. 37] it was shown that the requirement of multiplicative renormalisibility constrains 



the construction of possible transversal parts of the e.m. vertexfunction. 
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4-3 Pion Nucleon Form Factor 



The pion nucleon form factor, which enters into pure hadronic models usually as an input 
parameter, has been much debated in the last few years. Especially the value at the pion 
mass shell is not directly accessible in experiments. For a recent analysis obtained from 
different experiments (ttN scattering, NN scattering, ...) see ref. p9[. In this respect there 



is certainly the need to calculate this observable from a more fundamental quark level. 
Within our approach the pion nucleon form factor is obtained by coupling an external pion 
current to the diquark-quark bound state using again Mandelstam's formalism. Because 
of parity conservation one notes that there is no contribution from the loop integral 



(52 ) where the pion would couple to the + diquark. So the relevant nucleon current is 
completely carried by the quark part (|5l|) . In order to specify the vertex function for 
the pion-quark coupling we use spontaneously broken chiral symmetry and Goldstone's 
theorem ||40|| : In the chiral limit the vertex function, which is nothing else than the pion 
Bethe-Salpeter amplitude, is proportional to the scalar self energy of the quark. This is 
a consequence of the fact that in the chiral limit (vanishing quark current mass) where, 
provided the pion mass vanishes (Goldstone's theorem), the pion Bethe-Salpeter equation 
is equivalent to the quark Dyson-Schwinger equation. Therefore the proper choice is 

K{Pt = Q 2 = ^p 2 )=^l,^f±T\ (65) 

where B(p 2 ) is defined in eq. (^) and the experimental pion decay constant /V = 0.093 
GeV is used. Although this choice of the pion-quark vertex is valid only at the pion mass 
shell Q 2 = 0, we assume that the offshell amplitude does not vary too much with the pion 
momentum and therefore the most obvious generalization [41] 



C(Q V) = *7 5 y j-(B((p - \Q) 2 ) + B(( P + \Qf)) (66) 

of the onshell vertex function is allowed at least for a few hundert MeV around the onshell 
point. 

When inserting the pion quark vertex function into eq. (|5T| ) and evaluating the loop 
integral we determine the pion-nucleon form factor. Furthermore since the pion-nucleon 
matrix element is parametrized as 

UQ 2 ) = (NiPj^S^lUQ^lNiP^S,)) 

= A + (P / ,5 / )[ 7 5^VTv(g 2 )r a ]A + (P J ,^). (67) 

the extraction of g^NNiQ 2 ) is straightforward. 



4-4 Axial Form Factor 



Finally we calculate the axial form factor Qa{Q 2 ) of the nucleon. Due to parity, we again 
have to consider only the quark part of the nucleon current, since an axial current does 



22 



not couple to the scalar diquark. To specify the relevant vertex function, where the axial 
current interacts with the quark, we use the chiral Ward identity which reads 



Q,r 5 ,( P + Q, P ) = s-\ P + Q) 75 - i,s-\ P ). 



(6* 



This symmetry constraint can be satisfied with the vertex |40] 



' A(p 2 ) + A(k 2 ) 

% o 7m 



+ % 



(p + fc)„ 

p2 _ fc2 



(A(p 2 ) - A{k 2 )) 



7P + jk 



-i{B{p 2 ) + B{k 2 )) 



7 5 y, (69) 



that is totally determined by the quark self-energies. Although the axialvector vertex 
function looks, up to the isospin generators and the explicit 75, very similar to the Ball- 
Chiu vertex ([57]) (used in the case of a vector current coupling), there is one important 
difference: In the last term the scalar quark self-energies B(p 2 ) and B(k 2 ) add up, which 
in the limit Q = p — k — > leads to 



Qll TUp,k)^r a B(p 2 ) l5 . 



(70) 



When comparing ([7(]) with eq. (|65|) , the pion-quark vertex, one observes the exact Gold- 
berger Treiman relation for the quarks in the chiral limit: 



hm^ M r^(p,fc) = /^(0,p 2 ). 



(71) 



Stated in other words, in the chiral limit the axialvector quark coupling for vanishing 
momentum transfer is completely dominated by the pseudoscalar coupling to a massless 
pion (conserved axial current). 

To examine whether the phenomenologically much more important Goldberger- Treiman 
relation on the nucleon level 



M 

QttNN — 9A~jr- 



(72) 



is fulfilled, we have to determine Qa{Q 2 ) by taking the matrix element of the axial current 
between the diquark-quark bound states. After using the axialvector vertex function 
to evaluate (^) and comparing the result with the familiar decomposition 



J a M) = {N{P f ,S f )\fJQ 2 )\N{P^S)) 



A+(P f , S f )^g A (Q 2 ) + Q,g P (Q 2 )} l5 A + (P t , S,) 



(73) 



of the matrix element, the invariant form factors are immediately accessible. Note that 
besides g A {Q 2 ), the axial form factor of the nucleon, there exist another form factor 
gp(Q 2 ), the induced pseudoscalar form factor with a pole at the pion mass shell. 
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4-5 Discussion of Numerical Results 



In this subsection we finally show and discuss our results for the nucleon form factors. 
We will concentrate again on three different cases: Using type I propagators, the solution 
of the corresponding Bethe-Salpeter equation and the appropriate vertex functions; using 
type II propagators with d = 1 and with d = 10. The so obtained electromagnetic form 
factors of the nucleon are shown in fig. (^). 

F/(Q 2 ) F e N (Q 2 ) 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Q z /GeV 2 Q z /GeV 2 

Fig. 7. Electromagnetic form factors of the nucleon. The nucleon mass in all cases is chosen to 
be M = 1.9m q and the amplitudes entering the calculation of the matrix elements are obtained 
with A = 2m q . 



First of all, we mention (as can also be seen in the figures), that it is not possible in 
general to obtain exact charge conservation. This is due to the impulse approximation, 
where the coupling of the photon to the exchanged quark (and also to the diquark-quark 
form factor) is neglected. Whereas in a Salpeter approach charge conservation can be 
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shown analytically and it is well known that the impulse approximation is sufficient 

to guarantee charge conservation of mesons within the Bethe-Salpeter approach p2|, in 
our fully relativistic nucleon calculation the situation is more complicated (see also ref. 
29fl). But it can also be seen in the figures that the violation of charge conservation AQp^ 
is at the most only a few percent (due to isospin symmetry AQn = —2AQp). In this 
respect the contribution of the neglected exchange graph is more important for the case 
of a minimal screening of the mass poles (type II with d — 1). 

Furthermore the electric form factor of the proton in calculations of type I and type II 
with d = 10 are in reasonable correspondence with the empirical dipol fit 



F e P (Q' 



1 



Q 2 /0.71GeV) 



(74) 



The result obtained with d = 1 falls short in this respect: The variation of Ff(Q 2 ) with 
the photon momentum is too weak. This feature is also present for the other form factors. 
Obviously the modifications of the propagators with a pure exponentional factor severly 
changes the spacelike observable nucleon properties. 

The electric neutron form factor shows a reasonable course. Note that a sensible descrip- 
tion of neutron properties was one of the motivations for the diquark concept [ 43|| , having 
the following intuitive (nonrelativistic) picture in mind: A lighter quark circling around a 
heavier scalar diquark inevitably leads to a negative mean square radius of the neutron. 
It is now satisfying to observe that our relativistic approach, although working with equal 
masses of the constituents, confirmes this intuitive picture; it is not trivial at all, that 
the quark and diquark contributions give rise to a qualitatively correct behaviour of the 
electric neutron form factor. 



Nucleon form factors 




type I 


type II, d — 1 


type II, d = 10 


Exp. 


Qp 


1 +4.1 • 10~ 3 


1 + 2.85- 10~ 2 


1+4.7- lO" 3 


1 


Qn 


-8.2 ■ 10~ 3 


-5.70 ■ 10~ 2 


9.4- 10~ 3 





Hp 


1.32 


1.68 


1.34 


2.79 




-0.64 


-0.77 


-0.65 


-1.91 


9a 


1.39 


1.31 


1.41 


1.25 


9ttNN 


10.59 


14.17 


10.89 


13-14.5 



Table 3 

In this table we show the basic results for the form factors. Note that all magnetic moments are 
given in in units of e/2M. All displayed observables are obtained in calculations, where Si(p) 
and 5*2 (p) were included. The number of considered Gegenbauer polynomials was terminated 
after convergence was achieved. 



Going to the Salpeter limit we also get charge conservation numerically. 
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Fig. 8. The axial form factor gA(Q 2 ) of the nucleoli and the pion nucleon form factor g^NNiQ 2 
are shown for type I and type II with d = 1, 10 calculations respectively. 



Axial form factor g A (Q ) 

type II, d=10 




Pion nucleon form factor g„ NN (Q ) 

type II, d=10 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Q 2 /GeV 2 




Fig. 9. The axial form factor gA{Q 2 ) and the pion nucleon form factor g^NNiQ 2 ) are shown for 
a type II with d = 10 calculation, including only S\{p) or including 5*1 (p) and S^p). 

For both proton and neutron, the magnetic moments are too small. This is, however, not 
surprising because large contributions coming from 1 + diquarks are expected for these 



observables ( ||32|| , f44j|). Nevertheless, the calculated values have magnitudes expected 



from a calculation involving only + diquarks p0| , |31 . 



Before discussing the axial properties of the nucleon, we mention the following behaviour 
of the form factors: When A is decreased, leading to narrower amplitudes in momentum 
space (compare fig. (f|)), the electric form factors for the type I calculation become steeper 
and approach the emperical dipol fit. While such a behaviour is expected near threshold 
|f45|j we do not observe this feature for type II calculations. In these calculations the form 
factors (at leat at small Q 2 ) are very insensitive on A and the actual value of the bound 
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state mass. 



Our results for the axial and the pion-nucleon form factor are shown in fig. (^). Again, 
we observe for a type II with d = 1 calculation a relatively small variation of the form 
factors with Q 2 . Furthermore for the calculation including type II propagators with d = 10 
the similarities to the type I case can be seen. A striking feature, however, is the good 
agreement of qa = 9a{Q 2 = 0) and g n NN = 9-kNn{Q 2 = 0) with the experimental values. 
Qa tends to be slightly too large whereas the pion nucleon coupling constant obtained 
with type I and with type II (d=10) is a little bit too small. Nevertheless, in all cases 
a qualitatively right behaviour is obtained. This is astonishing, because also for these 
observables (in analogy to the magnetic moments) sizeable contributions from 1 + diquarks 
are expected. Note that in ref. J^J also a relatively large value of qa has been obtained, 
including only + diquarks. 

By comparing the calculated values of g n NN with the ones obtained form the Goldberger 
Treiman relation (72) we observe a violation of the Goldberger Treiman relation up to 30%. 



This failure also signals the necessity for the inclusion of 1 + diquarks in the calculation. It 
is furthermore to be seen how much the exchanged quark contributes to these observables. 

In figure @ we show the influence of the lower component S2 (p) of the nucleon Dirac 
field by comparing the results obtained with all amplitudes with the result obtained with 
Si(p) alone. It is interesting to note that the values of g n NN{Q 2 ) at small values of the 
momentum transfer are influenced by S^p); while the values of Qa{Q 2 ) at larger values 
are slightly more affected. Nevertheless we conclude that having a large lower component 
is probably only due to the representation chosen for the Dirac matrices. They do not 
have a significant effect on observables studied so far. 



5 Conclusions and Outlook 



In this paper we have developed a covariant diquark-quark model suitable for the calcu- 
lation of nucleon observables. As discussed, confinement is put in by a modification of 
the free (tree-level) quark and diquark propagators. The structure of the Bethe-Salpeter 
equation, describing nucleons as diquark-quark bound states interacting through quark 
exchange, has been elaborated using an appropriate decomposition of the Bethe-Salpeter 
vertex function in the Dirac algebra. Furthermore, the coefficient functions of the cor- 
responding Lorentz tensors have been expanded in hyperspherical harmonics, exploiting 
an approximate 0(4) symmetry (which would be exact if the exchanged particle were 
massless). After discussing the various numerical methods we presented our results for 
the calculation including + diquarks. The advantages of the confining propagators are 
clearly seen in the absence of unphysical thresholds. 

The numerically obtained nucleon vertex function has been fitted to a simple analytic 
form and then used to calculate nucleon matrix elements and form factors. In particular, 
we considered electromagnetic, the axial and the pionic form factors of the nucleon. While, 
due to the various simplifications and approximations, a quantitative agreement of the 
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results with experimental values could not be expected, we nevertheless got reasonable 
results. Furthermore, the large lower component of the nucleon Dirac field was shown to be 
merely an artefact of the chosen spinor representation. Although, the various observables 
displayed that the developed picture of a nucleon is still oversimplified, they nevertheless 
give us confidence that we are on the right track in obtaining a reasonable and trustable 
baryon model. The next step of improvement will be the inclusion of 1 + diquarks in the 
Bethe-Salpeter equation and also in the nucleon matrix elements. Furthermore we plan 
to determine the internal structure of the diquarks (described crudely by the parameter 
A) by a microscopic diquark model. 

In order to fully exploit the advantages and possibilities of our covariant and confining ap- 
proach, the model will soon be applied to other processes. In particular, we will investigate 
the different observables associated with the reactions p + 7 — > K + A and p + p —>■ K + A 
which are measured at ELSA ||6| and COSY [(47], respectively. 



While in most calculations of nucleon structure functions the distribution functions of the 
constituents are merely parametrized (see e.g. ||48||), our approach offers the possibility to 
determine, in a first step, these distributions, after the constituents are specified. In ref. 



24 1 this was investigated, assuming the scalar diquark to be a spectator and the photon 



interacts only with the quark. 

We conclude, that the reported studies are a good starting point for further investigations, 
which certainly will improve our qualitative and quantitative understanding of the baryon 
structure. 
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A Dirac Decomposition 



We choose the following ansatz for the amplitude x M 

X»(P,p) = A 1 (P,p) l5 P»ZA+ + A 2 (P,p) l5 P»A + 
+ S 1 (P,p) TB p£EA+ + B 2 (P,p) l5 p»A+ 

+ c 1 (p,p) l5 t(pn-Y)^ + + c 2 (Pph^(pn + ^)~A + (a.i) 

In the rest frame of the bound state, P = (0, 0, 0, iM), we obtain by using 



.0 0/' Kpa- 0/' 

the following expressions: 
(paA^P.p) 



28 



and 



X 



((p^pB^p) 
V pB 2 (P,p) 



+ 



crC^p) 
<T(p<r)C 2 (P,p) 



(A.4) 



Using the redefinitions 



Bi — Bi — Ci, Ci — Ci 



turns eq. ( |A.1| ) into eq. (pij). In the rest frame this leads to 



X 



p(pcr)B 1 (P,p) + i{a x p){pa)C x {P,p) 
pB 2 (P,p) + i(crxp)C 2 (P,p) 



(A.5) 



Thus in the rest frame the meaning of the amplitudes is therefore as follows: the am- 
plitudes Ai are the time components, the components Bi are parallel to the relative 
three-momentum of the constituents, and the components in the third line in eq. (|2"T| ) are 
orthogonal to the momentum and the spin. 



B Analytic fits of the amplitudes 

Fit Quality of S 10 (p) 



1e-04 



8e-05 



6e-05 



4e-05 



2e-05 



Oe+00 




Fig. B.l. (5i (p) 



numeric 



S\o(p)fit) is plotted for the type II calculation with M = 1.9m q . 



After solving the Bethe-Salpeter equation the amplitudes S w (p) ... S 22 {p) are given only 
at the grid points of the momentum mesh. For further use in the calculation of matrix 
elements this a very unconvenient feature. Therefore we fit rational funtions to them. To do 
this the value of every amplitude at each grid point is therefore treated as a "data" point 
with a certain standart deviation o. We find that the following set of rational functions 
are suitable to take the behaviour at small p as well as the behaviour at large p properly 
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Fit coefficients for type I 




$10 


S u 


$12 




5*21 


$22 


a 1 


0.9992 


-0.6954 


0.03028 


-0.4191 


0.4113 


-0.0127 


b l 


0.3107 


0.3028 


0.4617 


0.2682 


0.1937 


0.4090 


a 2 


-0.0002 


0.2929 


0.0134 


-0.1841 


-0.3939 


-0.0058 


b 2 


7.8284 


0.3028 


0.1834 


0.0792 


0.1966 


0.2085 


x 2 ° 


6.786 • 10~ 5 


4.025 • 10~ 5 


9.474 • 10~ 9 


1.132 • 10~ 5 


2.347- 10~ 8 


1.207- lO" 10 



Table B.l 

The fit coefficients for a type I calculation with M = m g and A = 2m q are displayed. 



Fit coefficients for type II, d — 1 






Su 


Sl2 


-S'20 


5*21 


'S'22 


a 1 


1.6562 


-0.9674 


0.4988 


-0.4168 


0.3898 


0.5042 


b l 


0.2437 


0.2556 


0.1985 


0.0924 


0.1564 


0.1562 


a 2 


-0.6574 


0.6342 


-0.4686 


0.0463 


-0.3528 


-0.5069 


b 2 


0.3838 


0.3293 


0.2046 


0.0924 


0.1613 


0.1558 


X 2 ° 


5.305 • 10~ 4 


1.340 • 10~ 4 


5.331 • lO" 7 


1.035 ■ 10~ 4 


5.263- lO" 8 


3.594- 10~ 9 



Table B.2 

The fit coefficients for a type II calculation with d = 1, M = m q and A = 2m q are displayed. 



into account : 



S " W = S(i + W 

2 CL^ p 2 

2 a fc Q 

S21 (p) = 2^ 7Z 7T. 777^ 

2 A; 3 

^ 22(p) = S(i+"L 2 ) 4 



B.l) 

B.2) 

B.3) 

B.4) 

B.5) 

B.6) 
B.7) 
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Note, that the momentum p is taken in units of the quark mass. The fit coefficients 
a^,b^] i,j = 1,2; k = 1,2 are then determined independently for each amplitude by 
performing a x 2 fit. To judge the quality of this procedure we show in fig. ( [B.l| ) the 
squared deviation 

^ 2 = (5*10 (P) numeric ~ 5*10 (p) fit Y ■ (B-8) 

Suming up the squared deviations then leads to 

kmax kmax 

° 2 X 2 = A 2 = ( s w(k) nurneric - S w (k) fit ) 2 , (B.9) 

k=l k=l 

i.e. the merit function x 2 times the squared accuracy of the "data points" . In the tables 
( |B.1| ) and ( |B.2| ) we show the fit coefficients which reproduce the plots for the type I and 



the type II with d = 1 amplitudes given in figure (0) within very small deviations. 



C Evaluation of the loop integrals 



In Section 4.1 we derived the two loop integrals building up the nucleon current. In eq. 
& 

m 2 ) = J ^m,Pf)s(p + K(p + ,p^s(p^D(p d )x(p t , Pl ), (c.i) 

the momenta are definded as 
P± = P + \Pi + \Q ± \Q 

Pd = -p+lp t (C.2) 
whereas in eq. 

KiQ 2 ) = J ^x(^.P/)^(P+)rJ(P+.P-)^(P-)^(P < r)x(fl,ft)- (C3) 
the momenta are chosen as 

P± = ~P + \Pi + \Q ± \Q 

p g = P+i]P t . ' ' (C.4) 
As stated in the text, to evaluate these integrals it is useful to work in the Breit-frame 
Q=(0,0,Q 3 ,0) 

= (o,o,-ig 3 ,V M2 + i^ 2 ) 

Pf = (0, 0, 1Q 3 , ijM* + lQ 2 ), (C.5) 
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and further to express the loop momentum p in 4-dimensional spherical coordinates 

p = p (sin0 sin 9 cosip, sin0 sin 9 simp, simp cos 8, cosip) . (C.6) 

While the <p integration is trivial, the integration over p, and ip has to be performed 
numerically. 
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